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diffusion. We observe a transition from monotone profile to non-monotone profile. This transition 
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I. INTRODUCTION 

The anomalous properties of the state with current in 
one-dimensional integrable system has attracted consid- 
erable interests. Especially, the heat conduction in one- 
dimensional systems is a long -standing problem. It has 
been expected that the existence of conserved quantities 
implies anomalous conductivity of the heat current Q . A 
large number of numerical investigations have been done 

It is natural to consider the state with current as 
the non-equilibrium steady state (NESS), i.e., the state 
which is asymptotically realized from the inhomogeneous 
initial state ||,[Q. In other words, the NESS is the state 
at the convergent point. There reported many works on 
the NESS itself. However, the relaxation process to the 
NESS is not much investigated. After a large but finite 
time, what kind of profile does the observable quantity 
show? In the previous paper @], stating from the inho- 
mogeneous initial state (the temperatures of the left and 
the right are different), we obtained the homogeneous 
NESS for the transverse AA-model. Investigating the 
profiles in the intermediate time, we would obtain the 
diffusion of the temperature profile. Furthermore, in the 
integrable system, the existence of the conserved quan- 
tities may have a significant influence on the relaxation 
process. 

In this paper, we investigate the profile at the large 
intermediate time using the transverse XA-model. The 
chain is initially divided to the left and the right, and 
kept at different temperatures. The problem of inter- 
mediate profile with inhomogeneous initial condition in 
the AA-model was first studied by Antal.et.al. in 
They considered the magnetization profile in zero tem- 
perature with the reversed external field, the situation 
which makes the calculation simple. It was shown that 
the magnetization profile shows the scaling property, 
m(x,t) pa &(x/t). They calculated the explicit form of 
the scaling function $, which has the flat part around 
the origin. We investigate this problem with the aid of 



the C* -algebraic argument. Modifying the argument of 
Ho and Araki , we develop a method which is applica- 
ble to the present situation. The C* -algebraic argument 
makes the calculation much simpler, and enables us to 
consider the more general situation, including the finite 
temperature case. 

We obtain the scaling limit m(x,t) ss $(x/i). For 
the situation of Antal.et.al. Q, we reproduce the same 
result. We further consider the finite temperature case. 
In the finite temperature case, we obtain a remarkable 
dependence of the magnetization profile on the strength 
of the external field and the temperature. When the 
external field is large, or the difference of the temperature 
is small, the profile varies monotonically. On the other 
hand, when the external field is small and the difference 
of the temperature is large, the profile is not monotone 
and has two extremum points. This feature is absent 
at the zero-temperature. That is, this phenomenon is 
due to purely thermal effect. This can be explained by 
the velocity distribution. In section [n],we represent the 
model and the initial condition. In section III, the sca ling 
property of correlation function is derived. In section IV, 
we investigate the interesting property which the profile 
reveals. 



II. THE MODEL 

The Hamiltonian we shall consider has the form, 

oo oo 

n = i E K< + i+ + l E <>(!) 

n— — oo n— — oo 

where cr" (a — x, y, z) is the a-component of the Pauli 
matrix at the site n. This Hamiltonian is called the trans- 
verse AA-model. The Hamiltonian is written by the 
fermion operators using the Jordan- Wigner transforma- 
tion il: 



H 



1 v-~>00 f . + 

'2 £n=-oo [ a n+l a n + <On+l 
iE~=-oo (H«n-1), 



(2) 
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where a n and a\ are the fermionic annihilation and cre- 
ation operators on the nth site. 

Owing to the bilinear form of the Hamiltonian, the 
dynamics of the many body system can be written by 



2 



the dynamics of the single particle. So, we only have to 
consider the dynamics of the linear combination of a\ ; 

oo oo 

t(/)= e /(0«l> E i/(oi 2 <^- (3) 



l= — oo 



L = — OD 



Here, I denotes the l-th site. The summable sequence 
{/(01 can be interpreted as the wave function of the sin- 
gle particle which is living on the lattice. They construct 
the one particle Hubert space, with the inner product 



(/is> = E /(Os(o- 



l — — oo 

As usual, we define the norm of / as 



\.f\\ = E i/(oi 2 



\Z— — oo 



In the one particle Hilbert space, the Fourier transfor- 
mation is defined as 

f(k) ee J2 f(n)e' mk , f(n) = i- f f(k)e~ mk . 

n— — oo ^ 

By the Hamiltonian (||), a^(/) evolves as 

at(/)-at(e <th /), (4) 

where e 4 *' 1 is the dynamics of the single particle, repre- 
sented in the Fourier representation as 

e^/(fc) = e - l * (cosfc - 7) /(fc)- 
Hence, in the coordinate representation, we have 

(e lth f){n) = — [ dke-^ cosk -^e mk f(k). (5) 

From now on, we use the following notations. Usually, 
the Heisenberg representation of the observable A is 

e UH Ae~ ltH . 

In stead of this, we use the notation at (A); 



UH Ae -itH 



a t (A). 



In this notation, the dynamics of the single particle [f| is 
written as 

a t (at(/))=at(e«V). 

We also use the unusual notation about the states. 
Usually, the state is given by a density matrix p, and 
the expectation value of A is given by 

TrpA. 



In stead of this, we use w; 

TrpA <-kj(A). 

These notations are introduced because the usual nota- 
tions are mathematically ill-defined in infinite system. 
However, there is no inconvenience in interpreting them 
in the usual sense. 

The initial state loq we consider is inhomogeneous. To 
define it, we divide the chain to the left and the right. 
The left (n < 0) side is in the equilibrium under mag- 
netic field 7_ with inverse temperature /?_. We denote 
the state by w_. On the other hand, the right (n > 1) 
side is in the equilibrium under magnetic field 7 + with 
inverse temperature j3+ . We denote the state by w+. The 
expectation value of a' (g)a(f) in cj_ , uj + are expressed as 

Lo_(a\g)a{f)) = - f dk P -(k)g-(k)J_(k), 



oj + (a){g)a(f)) = - f dkp + (k)g+(k)f + (k), (6) 



where p± is determined by the valuables (3± , and 7± as 

1 



p r (k) 



I -|- e -/3 r (cos(fc)-7 r ) ' 

We used the Fourier-sine transform; 
o 

f-(k) = -i /(n)sin((n-l)fc) fce[0,7r], 



2i 



dkf- (k) sin((n - l)fc) n < 0, 



/+(*;) = -» /(n)sin(nfe) fee[0,7r], 



n=l 



2i [* 

f(n) = — / dkf + (k)sm(nk) n > 1. 



Because of the quadratic form of the equilibrium states 
w_ and the expectation values are evaluated by use 
of the Wick product with the two-point function. The 
initial state has then the following product form; 



uj (A- ®A + )=lu- uj+ (A + ) 



(7) 



where A_ and A + are arbitrary operators of the left part 
(n < 0) and the right part (n > 1) of the chain, respec- 
tively. 



III. THE SCALING PROPERTY 

Now, we investigate the asymptotic profile of physical 
quantities. Let us consider X n , some physical value lo- 
calized in the neighborhood of site n. The expectation 
value of X n at the time t is X(n,t) = u>o(ott(X n )), with 
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the notation of the previous section. The scaling prop- 
erty means 



X(n,t)*<S> x (^) 



f 



a 



i.e., for large t, the expectation value of X at the site vt 
is almost $>x( v )- Fig.l shows the situation. Each figure 
is the snapshot of the profile of some X, at the time to, 
3t , 10t . Note that the leaned area has the width to, 3to, 
10to, respectively. It shows that the wave of the diffusion 
expand with the constant velocity l.They can be written 
as 



3t 



to 

10t 



X(n,t) = $ x (^) 



Here, &x(v) is the scaling function represented in Fig. 2. 

To show this scaling property, we have to prove the 
following convergence to the scaling function <&x ; 

<5>x (v) = lim X(vt,t) = lim u} (a t (X vt )). 

t — >oo t — *oo 

As ujo is determined by the Wick product of the two 
point function, we only have to derive the limit for X — 
aHg)a(f) case; 



10to 

Fig. 1 The picture of diffusion. The x-axes is the site. 
Each picture shows the snapshot of the profile of the 
local physical quantity X at the time to, 3to, 10to. 
The wave of the diffusion spread with the velocity 1. 

T t is the operator which is defined as 

(Ttf)(n) 

( g{n,t) for \n\ < t(l - S) 

Ath 



(e ,th /)(n) for \n\ >t(l-5) 



(8) 



where 



u v {a) (g)a(f)) = $ a t( 9 ) o (/)0) = hm u; (a t (a ii {S vt g)a(S vt f))) 

where we used the m-sift operator S m on the one particle 
Hilbert space; 

(S m f)(n) = f(n-m). 



1 {tt 



9{n,t) = -l-t^l-{n/ty 



f (_ S mr 1 (n/t)) e ~ insin 1 (n/t)-ity / l-(n/t) 2 +in/4 



For example, to derive the asymptotic profile of the mag- 
netization m(v), we calculate 



f 



m(v) = hm w ( a t ( -a% t 



+/(±7r + sin^ 1 (n/i))e m(±Tr+sin 1 («/*))+«V 1 -("/*) 2 -W4 

(+forn < 0, - forn > 0). 

This term corresponds to the contribution from the mo- 
mentum k(n, t) where the phase velocity 



lim u>o 



(a t {a\ t a vi 



<f){k) — — cos k + 7 + — 



= lim ujo (a t (a f (S vt r]o) a (S vt r]o))) - 
= uj v (a f (r?o)aM) - ~, 

with rjo , a wave function defined as 

/ 1 n = 
[0 n f= 0. 

Using the C*-algebraic method, we can make the argu- 
ment generic. As an advantage, we can treat the finite 
temperature case. The argument uses the result of Ho 
and Araki R|. They calculated the following decompo- 
sition of (||) with /:analitic function; for |n| < i(l — <5), 
with S > 0, we have 

(e M f)(n) = (T t f)(n) + (Af)(n,t). 



is stationary i.e. 



4>' (k) = sin k + - = 

=*> k(n, t) = - sin _1 (n/t),±7r + sin _1 (n/t) 
+ for n < 0, — for n > 0. 



(Af)(n, t) decay as 



\(Af)(n,t)\<^- 



(9) 



(10) 



f 



Fig. 2 The scaling function 4>x(v) corresponds to Fig.l. In this 
figure, the :r-axes is not sites, but the scaling factor v. 
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with C 5 , some constant which is independent of n. They 
also showed that the contribution to c lth f from \n\ > 
t(l-S); 

Sa = ^S, E l(e i4 V)(n)| 2 -0 (11) 

|n|>t(l-(5) 

goes to zero as <5 — * 0. (|To|), ( pi] ) ensures the following; 

e l "7 ~ T t /, (12) 

for large t. We can see from (||), this means that at the 
site n, there is only the particle at momentum k(n,t). 

Now, let us return to the current problem. First, we 
define projection operators on the one particle Hilbcrt 
space. The coordinate projections P~,P+ are 



71 > TO 

f(n) n < m 
f(n) n > to 
n < to 



> . 



There is the following relation between S m and P ; , 



c p± _ p± c 



The velocity projection P„ , P+ are defined in the Fourier 
representation as 



[Pv-f](k) 



kel v 
/>) kel c v7 



where 



I v = {k e ( — 7r, 7r];« < sinfc} . 



and 1^ is the complement of I v . By the definition of T t 
(S), we have for |n| < i(l - 5), 

(T t P±/)(n) = (P^ t T t /)(n). (14) 

Next we calculate the asymptotic value of 

at(o t (S ot /))=ot(e fth S rt /). 

We claim that 



lim \\P^S vt e lth f - P+S vt e lth f\\ = 

t— *oo 



We have to calculate 



(e" h 5^/)(n) = -!- / d/br l * (cosfe - 7 V (n -" t)fc />), 

27T /__ 



which we obtain by replacing n of (||), with n — vt. Note 
that S m commutes with the dynamics e lth . Let us define 
[g] t as 



(S vt T t g){n) \n-vt\<t(l-S) 
(S vt e lth g)(n) \n - vt\ > t(l - S) 



Using the relation Jl5|), (|lj), we have for |n — vt\ < i(l 

5), 



(n) = {S vt T t P+ g){n) 



(15) 



= (5, f Pr„ t T t5 )(n) = (P -S vt T t g)(n) = (P " [ 5 ] t )(n). 

(16) 

Let us fix some (5 > 0. Then, in the asymptotic limit 
t — > oo , we have from (|l^) , 

\\Pol9]t-[P v + 9}t\\ 2 

E KPo-^e^^H-^e^^P+gWI 2 

|n— ot|>t(l-*) 

+ E l^o-[3M™)-[P+5]*W| 2 
|n— ut|<t(l— 5) 



< 2 



E |e"V»)| 2 + E le^ffWI 3 

|n|>t(l-<5) |ra|>t(l-<5) 
^2(p 5 (.g)+P 5 (p+.g)) . 

We also have in i — > oo limit, 
l|5 0t e*" fl -[ff] t || a 

= E |(^e^g)(n)-(M 9 )(n)| 2 

|n-ut|<t(l-<5) 

= E \(e* th g)(n-vt)-(T t g)(n-vt)\ 
\n-vt\<t(l—S) 

= E l(e^5)(n)-(T t .9)(n)| 2 



(17) 



\n\<t(l-S) 



= E l(4?)M| 2 < E 



S\2 



\n\<t(l-S) 

Hence, we have 



M<*(l-<5) 



t 2 



(18) 



lim \\P^S vt e lth g-P+S vt e lth g\\ 

< lim \\PvS vt e m g-Po[g] t \\ 

t—>oo 

+ lim ||P -[ ff ] t - [P+g) t \\ + lim || [P+g] t - P+S vt e uh g\\ 

< 2(B 4 (ff) + B s (P+g)). 



We used (p7|), fll8[ ) and the commutativity of P+, 5„t, 
and e It?l . Note that the first term limt_>oo \\Pq S v te tth g — 
P+ S vt e tth g\\ is (5-independent. So, taking <5 — > 0, we have 

lim ||P -5 t , f c' it ' l . 9 - P+S„ t e ith ff || = 0. (19) 
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Similarly, we have 



lim \\P+S v te uh g-P-S vt c uh g\\ = 0. 

t— »oo 



(20) 



Next, we substitute the above result to the initial state 
to derive the two point function oj v {a) {g)a(f)). For the 
purpose, we rewrite the initial state with the projection 
operators P±. Note that 



' 22, fn sin ( n - l ) k 



n<0 



^i(n— l)k 



-i(n— l)k 



2/ 



-~ [P -f(-k)i 



— ik 



P -f(k)e l 



Hence we have 



Hence, the explicit form of the two point function of 

w v (at (g)a(f)) is 

1 



u v (a< (g)a(f)) = 



2- 



u<sin k 



2vr 



dkp-(k)f(k)g(k) 
dkp+(k)!(k)g(k). 



>sin k 



Note that w v (<v (g)a{f)) is translation invariant. 

The velocity of the particle with momentum k is sinfc. 
We see that only the particles at temperature with 
velocity (v < sinfc < 1) contribute to the state uj v . This 
represents the situation that on the inertial system which 
moves with velocity v, quasi particles with velocity less 
than v in the left part of the chain go to left infinity 
and do not appear in the correlation function lu v . This 
feature is also the case for the particles with temperature 
V/3+ 



w_(ot(ff)o(/)) 

1 

47T 



(21) 



dk 



P -f(-k)P -g(-k) + P -f(k)P -g(k) 



e^P -f(-k)P -g(k) - e-™P -f(k)P -g(-k) 



(22) 



We are now interested in o-i_(a^(g t )a(/t)), with gt = 
S vt e lth g, ft = S vt e ith f. We 'can see from '((22]), 
that uj-(a^ (gt)a(ft)) is written with PQS vt e lth g and 
Pt\ S v tC lth f '. Hence, we can apply the asymptotic form 

&■ 

Applying the formula, we have 



lim w-(a}{g t )a(ft)) 

t — >oo 



= lim — 

t— >oo 47T 
Q -2ife p+ 



dk 



2P+(k)f(k)g(k) 



: P v + (k)P v + (-k)e 2Mk f(k)g(-k) 



e^P+(fc)P+(-fc)c 



-2ivLk 



f(-k)g(k) p-(k) 



Here we used the fact that the Fourier representation of 
S m is 



S m f(k) = e-" n \f{k). 

By the Riemann-Lebesgue Theorem, the second and the 
third term vanishes and we finally obtain 

lim w_(at(ft)a(/ t )) = / dkp.(k)J{k)g{k) 



Similarly, we obtain 



1 



lim oj + (a^g t )a(f t )) = — / dkp + (k)f(k)g(k) 



IV. THE ASYMPTOTIC PROFILE 

With the two-point function, we can calculate the 
asymptotic profile of physical quantities. The magne- 
tization profile m(v) is 

m ( v ) = ir- I dkp + (k) 

^ 7r Jv>sin k 



1 

"2n 



dkp-(k) - i. 



/ v<sin k 

The magnetic current profile is denned by 
J M (v)= lim u,(a t (J^)), 

t — >oo 

where J^ 1 = S-{S^ +1 — S*Sf l+1 is the magnetic current 
at the site n. Similarly, J M (v) is calculated as 



'"<•> = h 

1 



dkp+{k) sin fc 



i;>sin k 



dkp- (fc) sin fc 

' t;<sin ft 

Let us consider the zero-temperature case: each side is 
in the ground state (/3 + = /3_ = oo) with magnetic field 
7+, 7_ respectively. We further assume 7 = 7+ = —7-. 
This is the situation that was considered in 0, and we 
confirm their results, 

f 



m(v) = —m(—v) = < 
where 7 = sin 7rmo 

J M ( V ) = J M (-V) : 



< V < COS7TTO() 

arccos(u) 
-m H cos7rm < v < 1 



1 < u 



r 1 1 . 

—7 = — sm7rmo < v < costtitiq 

7T 7T 

— \/ 1 — COS 7T7T10 < V < 1 

7T 

1 < v 
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In the zero-temperature case, regardless of the strength of 
the external fields 7+ , 7_ , we can show that the magne- 
tization profile m(v) is monotone. The situation is classi- 
fied by the following conditions: (A) the absolute value of 
the external field: (i)|7_|, |7+| < 1, (ii)|7-| < 1, |7+| > 1 
or | 7 +| < l,|7-l > 1, (iii)|7_|,|7+| > 1 and (B) 
the signs of the external field: (i)sign(7 + ) = sign(7_), 
(ii)sign(7 + ) = — sign(7_). Hence, we have 3x2 = 6 
cases. By the explicit calculation, we can show that for 
all the situation, the magnetization profile m(v) is mono- 
tone. Fig. 3 shows the magnetization profile for each case. 
As seen in the following, the finiteness of the tempera- 
ture destroies this monotonicity when the external field 
is small. 

To concentrate on the thermal inhomogeneity, let us 
consider the situation 7 = 7+ = 7_ 7^ 0. Due to the 
non-zero 7, the spins have finite magnetization up to the 
value of the temperature. The left side of Fig. 4 shows 
the profile of the magnetization for various values of the 
external field and temperature. Fig. (a) corresponds to 



the case 7 = — 0.5,(3- = 10, (3+ 



1. It shows the 



non-monotone profile. As the magnetization profile is 
monotone in zero temperature /3_ = (3 + = 00, this is 
considered to be a purely thermal property. When we 
increase the strength of 7, the non-monotonicity is lost. 
Fig.(b) shows the 7 = — 1 case with the same tempera- 
ture: (3- = 10, (3+ = 1. We can see the monotone pro- 
file. On the other hand, decrease of the difference of the 
temperature also destroies the non-monotonicity. Fig.(c) 
show the 7 = —0.5 case with the different temperature: 
(3- = 2, (3 + = 1. The difference emerges also in the mag- 
netic current. The right side of Fig. 4 shows the profile of 
the magnetic current. In the case (a), the current takes 
the maximum value at two points. They corresponds to 
the two extremum point of the magnetization profiles. 
On the other hand, in the monotone case, the current 
takes the maximum at the origin v = 0. 

The property of the profile (monotone/non- monotone) 
can be explained by the velocity distribution. Firstly, 
recall that the velocity of the particle at momentum k is 



-1 1 



(a) 





(b) 



0.4 

i 0.3 
0.2 




-1 1 



(c) 




Fig.4:The profile of the magnetization and the magnetic current 
at finite temperature for various external fields. The left graph 
shows the magnetization profile, and the right one shows the 
magnetic current profle. In the left graph, x-axes is the scaling 
factor v and the y-axes is the value of the magnetization m{v). 
In the right graph, rc-axes is the scaling factor v and the y-axes 
is the value of the magnetic current J M (v). (a) 17 = —0.5, (3- = 10, 
13+ = 1, (b): 7 = -1, p- = 10, 13+ = 1, (c): 7 = -0.5, f3- = 2, /?+ = 1. 

sink. To see the dependence of the velocity distribution, 
let us consider the derivatives of m(v) and J M (v) with 
v. For simplicity, we restrict ourselves to the case that 
both of p-(k), p+(k) are continuous, i.e, the both sides 
are initially in finite temperature. In this case, we can 
differentiate m(v) and J M (v) with v for — 1 < v < 1. We 
have 



dm(v) 
dv 



dJ M {v) 



1 



\p+(v) -P-(v)} , 



\p+{v) -P-{v)] , 



dv Vl - v 2 
where p_(resp.p + ) is the velocity distribution, 



Pr(v) 



1 

2^ 



1 



+ 



From these expressions, we can see that the difference 
p+{v) — p-{v) determines the monotone/non- monotone 
of m(v) and j(v): if, for example, there exists a point 
where p+(v) > p-(v) changes to p+{v) < p_(t>), the 
profile is non-monotone. Fig. 5 shows p+(v) for 

the above mentioned situation. We can see the crossing 
only for the case of non-monotone (a) . 



DISCUSSION 



We have investigated the profiles of the magnetization 
and the magnetic current, in the intermediate time to- 
wards the non-equilibrium steady state, using the trans- 

£ versed AA-model. We have found an interesting prop- 

Fig. 3:The magnetization profile at zero temperature for various , , , , c , , , , c 

I 1 c u r„ • ^1 i- r x mi • erty: depending on the strength ot the external holds 

externai hefds. 1 he a>axes is the scaling factor v. 1 he y- axes is J ^ 

the value of magnetization. All the case is classified by (A)x(B): 

(A) the absolute value of the external field: (i)| 7 _|, | 7 +| < 1, 
(ii)|7-| < l,|7+l > 1 or | 7+ | < 1, | 7 _| > 1, (iii)| 7 -|, > 1, 

(B) the signs of the external field: (i)sign( 7 +) = sign(7_), (ii) 
sign( 7+ ) = -sign( 7 _). (a):A-(i),B-(i), (b):A-(i),B-(ii), 
(c):A-(ii),B-(i), (d):A-(iii),B-(i), (e):A-(ii),B-(ii),(f):A-(iii),B-(ii). 
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and the values of initial temperature, the profile shows 
monotone/non-monotone property. This emerges as the 
result of the initial velocity distribution of the right and 
the left side. If there is a crossing between two distri- 
butions, the profile becomes non-monotone. This initial 
velocity dependence is due to a fact that the transverse 
AA-model preserves the one-particle mode. Each par- 
ticle runs to the infinity with its own velocity. In this 
sense, the integrability affects the diffusion profile in an 
essential way. 

The derivation of the asymptotic profile is carried out 
by showing the equation (|l9|),(20). These equations are 
due to (|l4|); the fact that for each site, there is only 
the particle with specific momentum. The specific mo- 
mentum is the momentum where the phase velocity is 
stational (^). In other words, if the dynamics of free 
Fermion is asymptotically dominated by the stationary 
point, we would have the same property as in this paper, 
even if the dispersion is not cosine. 
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Fig.5:The velocity distribution of the right and left sides. The 
:r-axes is the velocity v, the y-axes is the distribution p(v). 
The solid line shows p-(v), and the dash line, p + (v). 
(a): 7 = -0.5,/J- = 10, 0+ = 1, (b): 7 = = 10, f3+ = 1, 

(c): 7 = -0.5, ( 9_=2 I /3+ = l. 
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